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PHYSTCAL SULENCES 


ENERGY OF 1908 SIBERIAN METEORITE EXPLOSION CALCULATED 
Beijing SHENGXUE XUEBAO [ACTA ACUSTICA) in Chinese No 3, 198! pp 142-147 


‘Article by Hu Xinkang [5170 9515 1660) of the Acoustics Inetitute of the 
Chinese Academy of Sciences: “Calculation of the F cgy Released by the Great 
Meteorite Explosion From Infrasonic Data"*] 


[Text] This article proposes the use of infrasonic waves 
generated by a big explosion in the atmosphere to calcu- 
late the energy released by the explosion. We used this 
method to calculate the energy released by the explosion 
of the large meteor close to the ground in Central Siberia 
in the Soviot Union on 30 June 1908 to be (2,321.5) x 10°" 
erg, and compared this with other computed results. 





I. Introduction 


At 0016.2 hours on 30 June 1908 (Greenwich Mean Time), a large meteor exploded 
violently close to the ground in the Tunguska region of the great forest in 
Soviet Central Siberia. This incident shook the world at the time. Related 
scientific research is still being carried out today based on the information 
obtained at the time. In recent years, research has become more profound in the 
attempt to understand the nature of this incident. 


At the time of explosion of the falling meteor, strong infrasonic waves were pro- 
duced. These infrasonic waves were recorded by the microbarometrographs (Show 
and W.H. Dines, 1905) that had just been successfully developed and that were 
being used at several localities throughout the world. But because the mechan- 
ism of propagation of infrasound through the atmosphere at the time was poorly 
understood, the research results obtained from the data varied greatly (1-6). 


This article proposes a method of calculation of the energy released by the ex- 
plosion of the falling meteor based on the theory of infrasonic waveforms (7-8). 
The accuracy of the results is onefold, the true value aay be between one-hail 
and two times the calculated value. From the 1950's to the 1970's, many methods 
of estimation have been proposed but almost without exception, they can only 


— — — —— —— — 


*Thie article was received on 13 June 1980 








assure an accuracy of about 5 times (the true value gay be between one-half 
and five times the calculated value) (9-16), 


We lave uscd our method to caleulate the amount of energy released by the ex- 
plosion of the Largest falling meteor as yet known to mankind and compared it 
with the results obtained by various previous method, 

Il, Computational Formula 

According to the theory of infrasonic waveforms, the ‘afrason.c wave pressure 
produced on the ground surface at a horizontal distance & trom an explosion in 
the atmosphere is given by the following equation (7, 8): 


Pam SPOR, @, to, 8, 0) (1) 


where sound pressure P is the sum of the n simple normal waves of sound pres- 
sure, ad the simple normal wave of sound pressure is given by the following 
eyuation. 


= | D.F(e)de (2) 


o, = [aed] [ 2k" anindanies, 
2] [38 (a? Cag YP )~*) * (3) 


The meaning of the parameters in equations (1), (2) and (3) is as follows: 


where 


R is the horizontal distance between the origin of explosion and the observa- 
tion point, zo is the altitude of the explosion, z is the altitude of the 
observation point, Y), is the energy of explosion, k is the numbe: of waves, 
Oo is the density, M J. the acoustic channel factor, ® is the azimuth of the 
observation point, w is the angular frequency, 49 is the ratio factor #1. 


Le, Pe» te ate constants: their values are repspectively 1 km, 1.61 x 24.45 x 
1: Pa, and t. = 0.48 s. 


i is the observation time. 


ot a large explosion in the atmosphere, the main frequency of the pressure 
wave received at a distant locality (for example, the magnitude of the main 
angular frequency of the infrasonic wave column received at a locality several 
thousand kilometers from the origin of the explosion equivalent to 1000 kT of 
TNT is below 1/100), therefore we have: 


ww’ & Cage, )~? (4) 


a 
e 

















Substituting this into equation (1) we obtain: 


Pm Rs Vers S* Glo, am) (5) 


where >a Glas) in equation (5) is of course a complicated function. 


It is a function of the .'rtitude of explosion, frequency and acoustic channel 
factor. Theoretical and experimental research shows that for such an explosion 
far below the subsonic channel altitude, the altitude of explosion ‘as almost 
no effect upon the oscillating waveform or the amplitude. The subsonic channel 
altitude ie generally about 25 kilometers above the ground. Therefore, within 
the range of 0 to 3 kilometers sbove the ground surface, we can regard this as 
satisfactory. Also, many people (such as Scorer, Pekeris. Narkrider and Pierce) 
have shown in their research that regardless of whether the model of mltipie 
atmospheric strata, double atmospheric strata or single atmospheric strata is 
used, the foremo t parts of the column of infrasonic waves (the first three 
cycles) are basically the same. This means that the acoustic channel factor 
(which is atmospheric structure) affects the foremost parts of the infrasonic 
wave column very slightly, and therefore, if we only take the various parameters 
of the foremost parts of the wave column, then the acoustic factor M can be 
approximated as a constant (11). 


Obviously, under the above conditions, >. G.(«, 2) is only a function of angu- 


lar frequency w (or period T). According to the results of processing a large 
amount of experimental ¢.ca (17), we have finally obtained a simple formula to 
calculate the energy: 


Y=C+P+ Ts Rif? (6) 


where Y--energy released by the explosion, in unit erg; 
(=-empirical constant, with a value of 8 x 10/9; 
P--sound pressure of the column of infrasonic waves, in unit Pa; 
T--main cycle of the column of infrasonic saves, in unit s; 
R--distance between the point of explosic™. and the point of observa- 
tion, in unit Mn. 


According to the large amount of practice of many years, the accuracy of the 
results of calculation using this method can reach within onefold. 


III. Study of the Energy Released by the Big Meteoric Explosion of 1908 





Figure 1 shows the center of explosion and the recording stations at the vari- 
ous localities in the world. The longitude and latitude of the center of explo- 
sion is E 101° 57°, N 60° 54°. The author has organized the data and information 
of the time which have been gathered in Table 1. The table lists 10 receiving 
Stations; actually, the reception at the Greenwich station is a composite of 

six independent receiving stations. It can be seen from the actual data that 

the information gathered by England and by the Soviet Union itself are the most 
complete. Figures 2 and 3 are two typical infrasonic wave column diagrams. 


























Figure 1. Schematic Map of the Great Siberian Explosion Area and the [nfrasound 
Stations on 30 June 1908. 


*The base gap was taken from the world map of political regions in the world 
atlas published by the Map Publishing Agency in 1978. 
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Figure 2. The Record of Infrasound Produced by the Explosion as Received at 
the Slutsk Geophysical Observatory Near Leningrad on 30 June 1908. 
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Figure 3. The Record of Infrasound Produced by the Explosion as Received at 
Greenwich (England). 


Table 1. Detailed Data of the Infrasounds Accordixg to the Central Siberia 
Meteorite on 30 June 1908 
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(8) teave SR DARE TRS, 
Key : 
(1) Station name (5) Sound pressure 
(2) Latitude (6) Main cycles (s) 
(3) Longitude (7) Azimuth (deg) 
(4) Distance between the point of (8) Batavia is the original English 
impact and receiving station (km) name of Djakarta 


The simple normal waves in Figure 2 are basically superposed, while the simple 
normal waves in Figure 3 have basically been clearly separated. The frequency 
scattering characteristic can be clearly seen. 


Using the data of the 10 stations of the world listed in Table i and formula (6), 
we can calculate the energy released by the explosion of this falling meteor 
from the data of each station. The total of the results of these 10 stations 
gives the average vaiue and the standard deviation. The results are listed in 
Table 2. Finally, we conclude that ** energy released by the explosion of 

this falling meteor was (2.3+1.5) x 107" erg. 
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fable -. Calculated Energy of the Big Siberian Meteorite Explosion 
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Key: 
(1) Station name (9) Potsdam 
(2) Energy (erg) (10) Schneekoppe 
(3) Station name (11) Zagreb 
(4) Energy (erg) (12) Greenwich 
(5) Irkutsk (13) Batavia 
(6) Slutsk (14) Washington 
(7) Copenhagen (15) Average value of energy” 
(8) Berlin (lo) “The error is standard deviation 


From 1925 to the present, many methods have alre dy been used to study the 
energy released by the explosion of this falling meteor, the details of which 
are listed in the appendix. The estimated value of the energy is from 10'’ to 
10° * erg. The range reaches 7 degrees of magnitude. The estimated mass is 
from 13 tons to 500 tons, a range of over 5 degrees of magnitude. 


se should point out one thing in our calculations: for the main cycles about 
»,000 kilometers apart, we have used the values of the Greenwich station. The 
reason we can do this is that at such distances, further frequescy scattering 
of each simple normal wave becomes very weak. Of course, with more complete 
data, this is unnecessary. 


APPENDIX 


The values of the energy released by the explosion of this falling meteor 
derived (or estimated) by the various methods used from 1925 to 1980 are as 


follows: 


The method of using the audible threshold of audible sound produced a 
power of >> 10'” erg/s; 


i 
. 


Comparison with the big explosion of the Krakatoa volcano in 1883 gave an 
energy of ¢ 10?! erg; 


}. Comparisor with the strong winds of 50 meters per second gave a power of 
> 5x 10 erg/s; 








10, 


Li. 


12. 


Comparison with earthquakes gave an ene:gy of « 10°! erg; 


Estimates based on the height of the columns of fire and the intensity of 
light at the time of explosion gave a power of 2 x 10'? to 102" erg /3; 


Estimates based on the area of forests and trees destroyed by the shock 
waves of the exp'osion gave an energy of » 4.4 x 10°! erg; 


Estimates based on atmospheric fluctuations in the Soviet Union gave an 
energy of » 0.5 x 10°” erg; 


Estimates based on atmospheric fluctuations in England gave an energy of 
> 3.2 x 10°° erg; 


Comparison with the series of atmospheric nuclear tests conducted over 
Novaya Zemlya in the 1960's by the Soviet Union gave an energy of about 
10*" erg in magnitude; 


Estimates based on the depth and the dimension of the meteor crater gave 
a mass of about 13 tons; 


Estimates based on the greatest possible velocity of a meteor of the solar 
sys*em entering the earth's atmosphere gave a mass of at least 30,000 tons; 


The Ukranian Academy of Sciences of the USSR calculated the mass of the 
meteor to be over 5 million tons on the basis of the distribution of iso- 
topic microscopic granules and power of the explosion. 
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PHYSICAL SC LENCES 


OCEANIC SOUND SIGNAL TRANSMISSLON CHANNELS STUDIED 
Beijing SHENGXUE KUEBAO /ACTA ACUSTICA] in Chinese No 3, 1981 pp 152-161 


Article by Chen Cene [7115 16469] and Xu Junhua [1776 0193 5478] of the Acous- 
tics Institute of the Chinese Academy of Sciences: "Using the Method of Corre- 
lation Between Signal Pulses To Measure the Time-Variant Characteristics of 
Oceanic Sound Signal Transmission Channels"* } 


‘Text]} This article presents a method of correlation between 
signal pulses to study the random time-variant characteristics 
of oceanic sound signal transmission channels. It derives the 
expression of the correlation function between periodic pulses 
transmitted through a sound signal transmission channel and 
gives the results of measurement of the correlation coefficients 
between the pulse transmitted through the sound signal trans- 
mission channel and the pulse in the reflection channel of a 
fixed target, and it gives a preliminary analysis of the modei 
of the communications channel and the results of measurement. 


I. Correlation Between Pulses and the Model of the Communications Channel 


l. Derivation of the Correlation Function of Pulses of Time Variant Communica- 
tions Channels 


In recent years, a lot of research concerning underwater acoustic communications 
channels has pointed out (1-3) that because of the effects of scattering and re- 
fraction caused by undulations of surface waves of the ocean, unevenness and 
unsmoothness of the layers of the bottom of the sea and unevenness of the aqua- 
tic media of sea water, the actual oceanic sound signal transmission channel is 
not only multiple but it is also randomly variable in time and space. We can 
regard the underwater sound signal channel between the receiving point and the 
fixed transmission point as a linear random variable wave filter whose time- 
variant characteristics can be described systematically by a correlation func- 
tion: 


Ryfiw,, 1, Qo, t) = Bl HCia,, 1,)H* (iw, t,)) (1) 





*This article was received on 18 December 1979. 

















where H (.0, ©) te the traasfer function of the time-variant syeten, 


The input and output relationship (Pigure 1) of the linear system can be repre- 
sented by 


Yuj= _ ACime, ONO: = \" Alt, XC = edt (2) 
where 
— | - ‘= 
Ait, a) 5 |, Mltw, He dw (3) 


The expression of the correlation function of the system's output Yt) is 
Re (tay te) = ELV Cn)VC,)) 


mB \\ ACs AC mes te DN Cte = DC, — te dn, 


k, © represent the labels of the sample points at times t, and t,. in the random 
process Y(t). If X% and h are statistically independent, then 


Ry (ty, 4.) = \\ El AC ay, EL XC — eg X Ct, — 4.) dng, 


= \\ El Au, ty aC, t.)) : Rt — = , dugde, (4) 


where k. is the correlation function of the input signal. 








Figure 1. Schematic Diagram of a Time-Variant System 


y ef) L adit) 


— F 








Figure 2. Schematic Illustration of the Response of the Channel System to 
Periodic Pulses 
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if a determined periodically repetitive signal X(t + mT) is sent at the 
system's input end, we can refer to the input and output timetable to define 
the range 


Nal) @ XCmT + 2) (aT, mT + T,) 
YO YOeT ean et) [mT + wel ea + T,+T,) 
a=, 1, 2, coe N = l, 


where T == the width of te pulse of the input signal; 
T, == time expansion caused by the communications channel. 
= 
Y,() -= ¥(%) - \, AC my, POL — = ty day 
= 
VAG) VC) | ACs tN Cty = te ing 


We then delay Yp (t) which has been received fireat to find the correlation 
function between the pulses Y. (t) and 3 (t). Let 


t *t, #1’ 


In equation (7), ©’ is the delay time; for convenience, the apostrophe will be 


omitted, 


Ringe r)= = \ Yu) ¥(y > edt 
'. 
- T} \| AC my, ty Cue, tye + EXC = ma X(ye t= M dupdagdty 


If the syvetem is wide-sense stationary, and if we assume that the response 
function of the communications channel varies relatively slowly with t, in 
the interval (0,T), then we can obtain the following by approximation 


Ric) \\ t(my, pa Ca,, + £) E ( Xin — a Xl +t — —2 dugdn. 


- I AC, PpaACa,, b> COR Ct > ay mm w, lege, 


t represent the average value over the area of integration. We calculate 
the statistical average of *, " on both sides of equation (9) 


E\Rit)) = E\\ A agt AC ay ER > ty — ty lagu, 


(5) 


(6) 


(7) 


(8) 





(9) 








and exchange the order of & and the integration symbol and take into con- 


sideration the facet that »\ te tve determined signal 


R(t /?= | gi W(y, i)ACu,, i * FIRS e * a = He dtgda, 


(10) 


it can be seen that when there ia 4 periodic determined signal at the input 
end of the system, the correlation function Ry (+) of Y(t) at the oucput end 


of the system ie the special case of equation (4). 


Substituting equation (3) inte equation (10), we have 


Rin-= + J 338 


α. 


fecause when t <0, h t) = h (up, te) = 0, and the zero limit of the 


integral can be changed to = @, then substituting equation (1) into the above 


equation 
- ref | Riuk was aw, tr) . Rr +> a ue, pe rerentere 
' deyptu day e, 


wh ere A we 4 — Ww ‘ 
* 


Also because of integration 


\ Rt > wy mw, ere dayda, 


— RiCr eh rereeedadu, 


= X(e,)e"**"8( ty — oe) 


Substituting equation (12) into equation (11), we have 


Rc) : i{ Rik oy, Qa, TN (ay 0 ay — hy heed, 


‘ae 
«* 


- ~ Rul aw. 0, t) Xe)" deo 
‘iy . = 


where X (so) is the spectral density of the power of the input signal. 


(il) 


(i2) 


(13) 

















If a sine filler pulse of frequency wo, amplitude A and width T. is added at 
the input end of the system, ; 


X(a) a (AE2) [( teh —eaTan )' 4 (S86e + eT aa | 


(wo = we) T pn (o> we)T py: 


When T, is sufficiently wide, the (sin x)/x function can be substituted by 
the 4 benctden. thus 


X(@) ~ (47s) [BC o = we) + BC > w)) 


(14) 
Substituting equation (14) back into equation (13) we obtain 
Rye) be (Ae) |” RuCsns Os Few — om) + Oe + so) ide 
= ° 
gi (AT PR alone 0, ©) + commer (15) 


From the above analysis we see that changing different wo and actually measur- 
ing R,(t), the correlation function Ry (w , 0, t) of the system can be obtained 
from equation (15). When Ry (wo, 0, 1) cSrrespnd ing to X(w) is a slow variant 
function of w, Ry (w, 0, t) in equation (13) can be taken out from the integral. 


Ry) = α. 05 #) |" XC)e~de 


= + Rul, 0, ©) Re) (16) 


where wo is the central frequency of the signal. 


According to equation (7), we can write t in R, as t = (nm = k)T + ty and omit 
the periodic term (n - k)T; thus equation (16) becomes 


Re) = + Reo, 0, €) + RA) 
. (17) 


In general, regardless of the form of the signal used, from the actually mea- 
sured correlation functions between pulses, we can obtain the system function 
Ry (eo, 0, t) that describes the characteristics of the communications channe! 
from equation (17). 

















/, Some Explanations of the Characteriatics of the Oceanic Sound Communications 


Channe! 


Ase everyone knowr, the actual oceanic sound transmission channel is very complex; 
it te affected by random and uneven factors that have different physical mechan- 
ieme and thus the sound signals pro, agated along the communications channel 
follow multiple paths and are randomly variant in time and space. Therefore, 
the system funetion that describes the characteristics of the communications 
channe! is also a random function of time and space variations; for example, 
hehtt, t, ©), H = Hile, t€, ©). If h or H ie a noninterfering pure random 
process (without any even portion), then the method of processing signals by 
matching the wave filter and the correlation between pulses will be completely 
inetftective., Conversely, if h and H are completely determined functions, this 
means the communications channel is determined; for example, the single path 
communications channel without distortion, the matching wave filter can reach 
its own most ideal processing increment of 10 log (2W,T), (W,--front bandwidth 
of the receiver; T--length of signal). For the stable multiple communications 
channel, the correlation coefficient between pulses (compared to the large sig- 
nal noise) will approach 1. But the actual situation is between the two, and 
the corresponding output is very closely related to the ocean situation, the 
ecean zone and the actual measuring time. The above shows the system function 
of the communications channel should consist of two parts. They are: 


Completely interfering determined process--thie is the "stable part" of the 
communications channel composed of the waveguide effect, stable sound velocity 
gradient, aad stable oceanic flow. In addition, the determined process also in- 
cludes the time-variant interference of the system function of the communica- 
tions channel. 


Noninterfering random time-variant process--this is the noninterfering portion 
of the system function of the communications channel composed of random undula- 
tions on the ocean surface, turbulence, and random motion of unevenly scattered 
objects. Therefore, the system's response function can be written as 


Aa, 1) = J Eghhu) + SEM a, 0) + SE, Me, 0) (18) 


where VEn VE. v E,are the uniformly weighted coefficients of the response 
function determined by the actual environment of the communications channel. 


The first and second terms of the above equation respectively represent the 
stable portion of the communications channel and the interfering components of 
the time-variant response function of the communications channel. They are all 
completely interfering determined transformations of the communications channel. 
The third term represents the noninterfering random transformation of the conm- 
munications channel. Thus, the correlation function between pulses of equation 
(4) can be written a- 
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Rin. | BUC — + JEM (my, i> VE img, ta) OV Bohl ie) 


+ V E\Mtigs ty) + VEC tes te))) « Relty = tgs te = eg digdn, (19) 
whereas 
El) BC Byki) + VEC, te) V/ Bote un) + Beings ))) 
m Eighdl y dad tig) + VBE, Cad ey AC tg ty) hel tg ACs t)) 
+ BL EAC, aw, )) (20) 


When the process is wide-sense stationary, equation (20) becomes 


El « ) = Eph my al ng) + V/ Bg B Chel my Cu.) + hele AC, )) 
+ ERC my, #,, t) 


™ AC tay tty) > ERC my, wy, ©) (21) 


In the equation, A(uk, u,) represents the sum of the first two terms, Rj, is 

the correlation function of the response function of the communications chan- 
nel. Substituting equation (21) into equation (19), and omitting the period 
term in the quantity of the correlation function of the input signal, we have 


R(r)= \ AC uy, RCT, > ty — ty ign, + \{ ER uy, &., 6) 


RA t, > ty — ty lgde, 


Appropriate adjustment of the delay time t, enables us to cotain 


Rr) = \( AC uy, RCO dgde, + i ERC ay, te, RAO duyde, 


— B+ Rr) (22) 





By normalizing R(t) we can obtain the correlation coefficient p(t) between 
pulses. ; ; 


(rym bs rr) (23) 


The constant b represents the “stable part” of the communications channel, for 
a definite tranemitting and receiving position, it is unrelated to time. The 
term y.(t) is the contribution of the interfering part of the system function 
to the correlation coefficient between pulses. 














Our Measurements of the correlation coefficient o,(t) between pulses of the 
sound propagation channel and the reflection channel of a fixed target in 
shallow sea show that the so-called "stable part" does exist, because p(t) 
lessens as t increases and approaches a certain “stable value" b, Because 
of the random time-variant property of the communications channel, 0, (t) 
drops as the correlative time increases, different dropping rates reflect 
different propagation characteristics, and thus it also shows that the meth, id 
of correlation between pulses is an effective means to study the time-variant 
property of the communicationg channel. The following will introduce the 
method of measurement, give an results of actua! measurements and give a 
preliminary analysis of the results. 


[l. Method of Measuring the Function R, (t) » to) of Correlation Between Pulses 
and the Results of Measurements 


l. Method of Measurement 


We mentioned above that the system's correlation function Ry(w, t);, t2) that 
represents the time-variant characteristics of the communications channel can 
be directly obtained from the system's output correlation function R,(ti, t2). 
In the following, we will introduce the method of measuring the Ry (ti, t2) of 
the one-way transmission channel and the two-way echo channel. 


(1) The measurements of R, (ty » to) of one-way transmission are shown in 
Figure 3. 

















Figure 3. Block Diagram of Measuring R(t), to) of a One-way Propagation 


Channel 
Key: 
|. Signal generator 6. Ocean surface 
transmitter 7. Inhomogeneous water cluster 
; Four-bit quantized compressed 8. Transmitting point 
time-correlator 9. Interior wave 
4. Time gate 10. Sea floor 
5. Recorder ll. Receiving point 
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The transmitter is located on the tranemission vessel, The transmission con- 
verter is placed at the botton of the sea at a fixed location, For reception, 
an underwater base fixed at the bottom of the sea is used, The tranemitter is 
controlled by the signal generator. <A frequency modulated pulse FM 100 and a 
32 x 20 coded phase modulated pulse are tranemitted at intervals of T = 5 
seconds. The pulse width is 640 ma. The pulsea signal travels via the under- 
water sound propagation channel to the 4-bit quantized linear time compressed 
correlator, and at the same time enters the pulse passage selection controller; 
the controller produces controlled square waves and groups the pulses received 
in order, every N number of pulses form a group, n= 1, 2, ..., N. The corre- 
lator stores the first pulse received as a reference signal aod correlates it 
with the second, third, ..., Nth pulse received. In this way, we obtsin Ry gir 
(t,), to), Ryqo(t), Cady wees Ryin (ti » ty) b= 1, 2, «ee, &k CL is ordering), 
to ~ ty = T, ty = ty @ 2T, woey Cy = ty = CN = 1)T, 


To eliminate the effect of amplitude undulations of the pulse. the correlation 
functions measured are converted to correlation coefficientr, i.e., 


Pyil tis Hr) ™ Ryal ty &)/[ Yala)¥ K4)) (24) 
where Yun) = Rin, = V Ryl0) 


Yin) =v Ry) = Vv Ryh0) 


The yi! + Pyi 2s see Oya (Nad) of each group are averaged according to their own 
quantity i = 1, 2, «.., k 


> 


bn = (>) Prin nmi, 2,*++,(N—1) (25) 


‘ 


to obtain the average value and the standard deviation for graphic represent a- 
tion. 


(2) The measurement of Ry (th, to) of the two-way echo wave is shown in Figure 4. 


The method of measurement is basically the same as that for one-way transmission. 
The difference is that transmission and reception are both fixed at the bottom of 
the sea at the same underwater base. The echo comes from a fixed target at a 
definite distance in a definite position selected. When a fixed target at a 
distance of 16 cable lengths is selected, pulses are transmitted at intervals 

of T = 10 seconds; when a fixed target at a distance of 56 cable lengths is 
selected, T = 20 seconds. 
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Block Diagram of Measuring Ry (th, to) of a Two-Way Propagation 
Channel 


6. Interior wave 
7. Ocean surface 

Four-bit quantized compressed-time 8. Inhomogeneous water cluster 
9 
0 


Transmitter 
2. Signal generator 
, 


. Sea floor 
. Transmitting and receiv ng point 


correlator 
4. Time gate l 
5. Recorder 





2. Results and Analysis 


Figures 5-7 show the curve of the correlation coefficients obtained by conver- 
sion from the actually measured R,(t;, t)). The actual conditions of measure- 
ment are labeled in Figure 5. "0" is the measured value, the solid line is the 
curve fitted by the least square method. 


Statistical analysis of the data shows the correlation coefficients between the 
pulses of the transmission channel and the fixed target reflection channel 
basically attenuate exponentially as time increases. Also, the time-variant 
characteristic of the communications channel can be explained in general by the 
following illustrative Figure 8. 


(1) For the determined ideal communications channel without noise interference, 
the correlation coefficients between pulses are unrelated to time, and 

o¥(t) = 1. 

y 


(2) For the determined communications channel (such as the stable and multiple 
path communications channel) with unrelated interference, the fixed correla: ive 
loss is produced by the noncorrelation of noise and mixed responses mixed in 
with the receiving signals; o,(t) does not change with time either, and sur- 
passes a certain constant (a J b). 


(3) For random time-variant communications channels possessing partial corre- 
lative characteristics, statistical analysis of a massive amount of data shows 
that the correlation coefficients between pulses drop exponentially as time in- 
creases; refer to equation (23). 


-aT 
o. ** ae +b 


























(s) (b) 


Figure 5. Time Variance of the Correlation Coefficients Between Signal Pulses 
in the Transmission Channel 


(a) The form of the signal of the 79.4.27 shallow sea propagation experiment, 
the transmitting and receiving points were fixed at a distance of 4 nautical 
miles: FM 100, t = 640 ms. 


§ The measuring point is the arithmetic mean of 22 sets of data, and the 
standard deviation is labeled. 


- The solid line is the fitted curve: o (t) = ae i+b a= 0.37; 
y a = 0,011; b = 0.6 


(b) The form of the signal of the 79.4.28 shallow sea propagation experiment, 
the transmitting and receiving points were fixed at a distance of 7 nautical 
miles: FM 100, T = 640 ms. 


© The measuring point is the average value of 6 sets of data (standard devia- 
tion oO. < 0.15). 


- The solid line is the fitted curve: o i) — i+ba. 0.436; 
y a = 0.016; b= 0.36. 


In the equation, a, a, b are constants and a, b are equivalent to Ey and E; in 
equation (18). E» represents the scattered part of the communications channel 
satisfying the wide-sense stationarity uncorrelated scattering (WSSUS) condition 
(1). They are determined by the actual environmental conditions. 





19 








e 
» * 
2* 
€ | ; i 
) * @ 





— > oa | 


r(s) 


Figure 6. Variation of the Correlation Coefficients Becween Pulses in the 
Two-way Communications Channel Along With Time 


The 79.5.7 shallow sea fixed target reflection experiment FM 100, 640 ms 





1. The target is 1.6 nautical miles from the base: The measured point "0" is 
the average value of 36 sets of data, the standard deviation is (0.13-0.18), the 
fitted curve is ° ae~"" +b 


a= 0.6; a= 0.0104; b= 0.36 


2. The target is 5.6 nautical miles from the base: The measured point "0" is 
the average value of 21 sets of data, the standard deviation is (0.13-0.21), the 
fitted curve is py (tr) =ae "*' +b 


a = 0.36; a = 0.0141; b= 0.30 


(4) It can be seen from the figures that the actual sound transmission channel 
also has the following outstanding characteristics: 


i. The first point of the actually measured correlation coefficient (i.e., the 
measured point corresponding to T seconds on the horizontal axis) drops by a 
larger scale, the latter neighboring points drop less in scale. We believe this 
is caused by the effects of random weighting of the signal by the communications 
channel and uncorrelated noise and mixed responses contained in the reference 
signal stored in the receiver. 


ii. It is also possible that a node Q exists. It can be derived from the 

first small point described in sections (2) and (4) that the time tg correspond- 
ing to Q should be the correlated radius of the random interfering background 
and it is determined by the bandwidth of the interference. 


iii. Results of measurements show that o,(t) will approach a stable va'ue b as 
time increases. This shows that regardless of how complex the transmission 








characteristics of the communications channel are, there will always be a 
stable even part. R. Thiele (3) measured the scattering function 6 full days 
apart at two fixed points (at a distance of 1.6 nautical miles apart) on 
Fehmarn Island in the Baltic Sea, and the similarity of the results also ex- 
plained this phenomenon, 
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Figure 7. Variation of the Correlation Coefficients of Pulses in the Two-way 
Communications Channel Along With Time 


The 79.5.11 shallow sea fixed target reflection experiment 
Target was 1.6 nautical miles from the base, FM 100, 640 ms. 


The measured point is the average value of 20 sets of data, and the 
standard deviation is labeled. 


- The solid line is the fitted curve. 
at 








p(t) = ae +b a=0.34 a#0,.021 b = 0.40 
| * Key: 
10 (2) (1) Correlation coefficient between pulses 
. al Oe PU TT (2) Noninterferi art 
(1) ng Pp 
e-/ 2S me (7) (3) Interfering part 


(4) Stable part 

(5) Random part 

(6) Correlative loss produced by 
noninterference of the interfering 
background 

(7) Increment of correlative loss caused 
by the time-variant characteristic 

nis of the communications channel 
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Figure 8. Illustrative Diagram of the Model of the Communications Channel 


T is the interval of consecutive pulses; "o" is the average value of the 
measured points; the solid line is the curve fitted by the least square method, 
y = ae +b. 














ill. Conelustion 


|. Actual Measurements show the oceanic sound transmission channel should con- 
sist of two parts, a determined transformation and a random transformation. 

The interfering part of the system function of the communications channel varies 
according to the differences in the actual ocean environment. 


2. The method of correlation between pulses is a good method to study the ran- 
dom time-variant characteristics of underwater sound communications channels. 

ine degree of the effect of the random time-variant part upon the transmission 

of sound in the communications channel can be expressed by the rate of exponen- 
tial drop of o,(t) with time. The variations of the interfering and noninterfer- 
ing parts of the communications channel can be directly obtained from measurements 
of the correlation coefficients. 


3. The system function Ry(t), t2, w) of the communications channel can be ob- 
tained from the measurement of the coefficients of correlation between known 
signal pulses. This method can also be used directly to increase the reliabil- 
ity of underwater communications. 


The above study of the oceanic sound transmission channel is preliminary and the 
data are insufficient; more profound work is still required. 


The work benefited from a great deal of discussion with Comrades Huang Zengyang 
(7806 2582 2543] and Zhu Ye [2612 6851]. Comrades Sun Zeng [1327 1073], Sun 
Puan [1327 4395 1344], Wei Xuehuan [7614 1331 3883], Zhou Wei [07.9 0251], and 
Zhou Guiqin [0719 2710 3830] participated in the experimental work. 
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SOUND WAVE PROPAGATION IN FAOLINITE SUSPENSION STUDIED 
Beijing SHENGXUE XUEBAO [ACTA ACUSTICA] in Chinese No 3, 198] pp 181-188 


[Article by Tang Yingwu [0781 2019 0710] of the Acoustics Institute of the 
Chinese Academy of Sciences: "Propagation of Sound Waves in Kaolinite Suspen- 
stons“*] 


Text) This article takes into consideration the effects 

cf che mutual reaction of fluid dynamics and the mutual 
reaction of electrical forces in an acoustic field of sus- 
pended particies, revises the equivalent dynamic density 
equation and the motion equation of a suspension, and derives 
a more perfect expression of the effective density of a sus- 
pension. The consistency between the results and the exper- 
imental data is good. 


l. Introduction 


The question of sound propagation in suspensions has been the subject of study 
by many acoustics workers (1) because of its importance in application and its 
difficulty of treatment by mathematics (when using sound scattering). it is a 
question that has not yet been appropriately solved in sound propagation theory 
during the past 60 years. Past theoret{._al work mostly studied the velocity of 
sound and the absorption coefficient in suspensions separately; only a few arti- 
cles have treated them together (2-4) (i.e., studying the question of sound 
propagation in suspensions), and the results obtained have neglected the mutual 
reactions of the suspended particles in the acoustic field. Such results are 
only appropriate for suspensions with a sufficientiy low volume concentration. 
Experiments prove that when the volume concentration is greater than 0.8 per- 
cent (5), the mutual reactions of suspended particles in the acoustic field 
will emerge. Experiments by Urick (6) and Hampton (7) showed that the mutual 
reaction of suspended particles in an acoustic field is an important factor 
affecting sound propagation in a suspension. Therefore, in discussing sound 
propagation in a suspension, consideration must be given to the mutual reac- 
tions ot suspended particles in the acoustic field. 





*This article was received on 27 July 1979. 
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[he Mutual reactions of suspended particles in an acoustic field mainly include 
iwo types: the mutual reaction of fluid dynamics (5) and the mutual reaction 
of electrical forees (9) (for pure kaolinite particles). This article will 
omit the effecte of the attenuation of .cattering and the attenuation of ther- 
mal conductivity upon the propagation of sound, and will discuss the question 
‘ff sound propagation in suspensions by considering the two types of mutual 
reactions described on the basis of Ament's (3) theory. The results obtained 
are compared with the experimental results, The theoretical and experimental 
resulte coincide very well, 


[l., Amalysis of the Mutual Reactions of Suspended Particles in an Acoustic 
Field 


Let us assume that the suspension consists of some electrically charged spheru- 
litic particles of pure kaolinite of even size suspended in pure water (abbre- 
viated water). When an acoustic field does not exist, the particles in the water 
are arranged in an equidistant cubic dot matrix. When there is an acoustic field, 
the wave line of the sound wave (the direction of propagation of the sound wave) 
is perpendicular to one side of thie cubic dot matrix. To facilitate analysis, 

“e number the particles and shade the plane of particles perpendicular to the 
wave line. The entire design is illustrated in Figure 1. From such an arrange- 
ment of particles, we can see the following: 





|, When there is an acoustic field, the suspended particles on the same shaded 
plane are relatively stationary; therefore electrical forces do not exist between 
these particles. 


’. Im an acoustic field, the mutual reaction of fluid dynamics existing between 
one suspended particle and the neighboring one on the same shaded plane is the 
major reaction compared to the mutual reaction of fluid dynamics between these 
particles and other particles. This is because for each particle, the greatest 
interference of “water flow” surrounding it is exerted by the several neighbor- 
ing particles. Because in an acoustic field, the “point mass" of water vibrates 
along the direction of sound propagation, thus when the wavelength of the sound 
wave in the suspension is very long compared to the linear dimension of the par- 
ticle; the “water flow streamline” between particles on the same shaded plane 
contracts iqward; this means that these particles exert the strongest resistance 
and friction on the flow of water, and thus, the mutual reaction of fluid dyna- 
mics between each particle and its neighboring particles is the strongest between 
those particles on the same shaded plane. Therefore, in this article we wiil only 
consider the mutual reaction of fluid dynamics between these particles. 


\. Im an acoustic field, the mutual reaction of electrical forces exists between 
the particle in the {th shaded plane and the several neighboring particles in 

the (j = 1)th and (j - 1)th shaded planes. In the following, we will analyze 

the electrical forces exerted on particle 0. When an acoustic field is not pres- 
ent, the distance between particle 0 and its closest neighboring particles 1 and 
1" is d' the distance between particle 0 and its next clesest neighboring parti- 
les 2, 4, 6, 8, 2", 4', 6' and 8 is v2 d; the distance between particle 0 and 
the neighboring particles further away, 3, 5, 7, 9, 3°, 5°, 7° and 9° is 3 4. 


Pr 








if we denote the electrical force exerting on particle 0 in the acoustic field 
defined by the particles |l-1' as f,, then from Pigure 1 we know that the elec- 
trical forces exerted on particle 0 in the acoustic field defined by particles 
d=2', &@4', and Be8' are directly proportional to {,, and the electrical 
forces exerted on particle 0 in the acoustic field defined by 3-3', 5<5', 7-7' 
and 9-9' are aleo directly proportional to f,, only the proportional coefficients 
are different. Therefore, in the acoustic field, the electrical foree Fe on 
particle 0 is 


PF eli em * ml, (1) 


Here, m, and m are two proportional coefficients to be determined, 


The electrical force upon the other particles in the acoustic field is also Fee 





Figure 1. Suspended Particles Are Arranged in a Cubic Dot Matrix in a Station- 
ary Body of Water 


lil. Effective Dersity 


Ament used the fol’ owing three equations (3) when he derived the effective 
density of a diluted suspension: 





em (l — Vey + Vea, (2) 
pe = (1 — Vow, + V om, (3) 
(4) 3 er cor’ pees — oye, ) = 04/3 er "CS — oF lye! o, — &)) (4) 
where 
— 
2 A 


(= 


Woliny |! * Torey) 
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» is the effective density of the suspension, v is the effective vibration 
velocity, Vol= & © my'N/)) te the volume concentration of the suspended parti- 
cle, y ie the radius of the particle, w is the angular frequency, » = n/6), 

o, ie the density of water, po» is density of the particle, » is the vise- 
eosity coefficlent of water, 1 = ¢=1, v) and vp represent the average vibra- 
tion velocity of the water and particle in the suspension, N is the number of 
suspended particles in a unit volume. The time factor of the acoustic field 
is chosen as At (t ie time), and is omitted, 


Based on the arrangement of particles described in section two and the mutual 
reaction among the particles in the acoustic field, we let 


ked << 1 (5) 


thus, afong equations (2) — (4), only equation (2) suite the present situation, 
and equations (3) and (4) must be revised before they can suit the present situ- 
ation. This is because equation (2) is a continuous condition, and regardless 
of whether or not, it is always correct. Equation (3) is the equivalent dyna- 
mic density equation, and equation (4) is the equation of motion of suspended 
particles in the acoustic field. They are only correct when gutual reactions 
among suspended particles in an acoustic field do not exist. When mutual reac- 
tion among suspended particles in an acoustic field exists, they are not correct 
any more. In the following, we will consider the mutual reaction of suspended 
particles in an acoustic field to revise equations (3) and (4). The ko above is 
the phase shift constant (the real part of the propagation constant) of the sus- 
pension, d is the distance between the centers of two particles. 





Firet we will discuss the revision of equation (3). 


The equivalent dynamic density ov in the suspension equals the sum of the motion 
(p" dy of the particles in a unit volume and the motion of water (ov), is the 
unit volume, i.e. 


ov = (ov), + (ov), (6) 

where 
(ov), = (1 — Vay — P., (7) 
(or), = Vee, = Vv \ F,- di + P., (8) 


Here, P.. is the viscous motion density, the density causing the loss in the 
motion of water due to viscous drag between the water and particles; y "a. 
is the motion density of electrical force, the density causing the loss in the 
motion of the particles due to repulsive forces in the acoustic field. For 
the arrangement of particles in Figure ', when equation (5) is valid, McCann 


(9) has already found 


= —Sd'iie, = > arp” * (9) 








Therefore Fe in equation (1) can be expressed as 


‘i= - ofall + om +> mm PV? 5 (10) 


where P is a constant, vo ie the inetantaneous vibration velocity of the 
particle, 


Substituting equation (10) into equation (8) and when t = ¢), v9 = 0, when 
t = to, Vo = Vo, thus 


Cow), = awe Lo — VY") > P, (11) 


where » Po/oy, H= CL + my + m)P. 


Substituting equations (7) and (11) into equation (6), we obtain the equiva- 
lent dynamic density equation in the suspension: 


pe = (1 — Vom + Valo — AVY de, (12) 


Now we will discuss the revision of equation (4). 


Because of the presence of mutual reactions of fluid dynamics between suspended 
particles in the acoustic field, the vibration velocity of the particles rela- 
tive to water (v> - v;) will increase by » times. Thus, equation (4) should be 
revised as 


© seri 7s — eye) = ar S — it ordeal vs = 1) (13) 


Here, » is a function of V., called the mutual reaction function of fluid dyna- 
mice between suspended particles in an acoustics field. If only the mutual re- 
action of fluid dynamics occurring between a particle and its neighbor on the 
same shaded plane is taken into consideration, we have (8) 


3 /6V,\" 1 (/6V,\|"* 
om [1-2 (S) ~i(&s)) (14) 


Also, because of the mutual reaction of electrical forces among suspended par- 
ticles in the acoustic field, the contribution of electrical forces must also 
be added to the left side of equation (13). Thus, equation (4) is finally re- 
vised to: 


‘(Co — HV )e; — 0) = (S — ie io; = 4). (5) 


Simultaneously solving equations (2), (12) and (15) yields the p expression for 
the effective density of the suspension: 


iofe — 1 — HY V1 — V,) 





p=, (16) 


bS —1(or +o — HV — Vie—1 — HVY))* 
whe re o.= (l= Ve + Vie. . 








Of course, when Ho = O (or P = 0), equation (16) is reduced to the expression 
(19) in Reference (10), When » = | and H = 0, equation (16) reduces to the 
resulte of Ament, or degenerates to equation (27) in Reference (11). 


LV. Propagation Constant 

hecause we do not consider the effects of sound attenuation caused by heat con- 
duction upon the propagation of sound in a suspension, the coefficient K of 
absolute heat compression in the suspension takes real values (4), and its 


expression is 


K= (1 - v) K) + Ve Ko, 


here Ky and Ko are the coefficients of absolute heat compression of the water 
and the pu e kaolinite particle. Substituting equation (16) into the expression 


t= w( Kp)” (17) 


vielde the sound propagation constant k in the suspension. 
In our case, the quantity ale —!i—AVi!yYVLI—V,) is much smaller than 1; 
therefore the real part Re(o) of o is much larger than p the imaginary part 
Im(o) of o, d.e., Reo) >> Im(o); thus equation (17) can be simplified to 

h = w(K Re(p))'*(1 + 1 Im(p) 2 Re(p)) 

- he > i, (18) 


From “9 we can obtain the sound phase velocity cp in the suspension: 


ce = 1 (KK Re(p))'*: (19) 


. is the sound amplitude absorption coefficient in the suspension; its 
expression is 


sg ep KA OW LL = V Xo = 1 = HV Pgs 
Re( pe) (oS + , 





(20) 


where 


Re (pop) =p = VAl — —— —— 





(oS) + 9° 


vu — ors c— AV e — |— Hv**),. 








Vv. The Expressions of the Sound Phase Velocity and the Abaorption Coefficient 


With equations (19) and (20), it is possible to theoretically calculate the 
sound phase velocity and absorption coefficient in a suspension, but it is not 
possible in the actual situation because the constant H is not yet determined, 
In the following, we have utilized the experimental data from Urick (6) to 
determine the value of the constant H. When the frequency is 1 megahertz, 
Urick measured the dependent relationship of the sound absorption coefficient 
upon the particles’ volume concentration in the pure kaolinite suspension; the 
measured results are shown in Figure 2, The physical parameters of the suspen- 
sion are: p); = 1.00 g/em’, oo = 2.65 g/em’, r = 4, 70 x 10~°em, wu = 0,01 /s 
equivalent to 20°C in temperature), cp 1, $2 x 10° cm/s, w = 27 x 10°, We sele 
selected various values for H and together with the above parameters we substi- 
tuted them into expression (20) for numerical calculation. We discovered in 
the results that the a-V. characteristic curve (the solid line in Figure 2) 
calculated when H = 0.7 and the experimental data coincide well (except beyond 
the range Ve » 0.35). In this way, we determined that H = 0.7. At this time, 
the expressions (19) and (20) of the sound phase velocity and absorption coef- 
ficient in the suspension can be written as 


C.- Kp _ KVLI—V Xo—1—0,7V! Yi br eo—0.7VM—V —— 2) | ae (21) 
a WS + (OS +o — 0.70% — Via — 1 = 0.7V™)}? 


K \'" oVLli—VXao— i —07V!)Se 


=e —— VS > (or — Vie =—1 =—0.7V™)) 








(22) 


Ordinarily, (rf +o¢—0.7V)) > VAe—1—0.7V)"), thus expressions (21) and (22) 
can be simplified to 








| . AVAL = V ho — 1 — 07 or +e — 0.7V3")\-" 
C, ™ (Ke, Ss + (or +o — 0.70) (23) 
K \* oVLl = V Xo = 1 — 0.70" ys 
ae ; a (OSY + (or eo —O.7VhY , (24) 


Expression (24) is the same as formula (16) in Reference (8). 


When the volume concentration of the particle in the suspension is very low 
(i.e., when V. << 1), expressions (21) and (2. can be further simplified to 





C, = (Ke, @ Cy (25) 
lo = —Y, ig = 1S 
c S+la+r) 
(26) 


Expression (25) is the famous Wood formula, while equation (26) is the same as 
the Lamb-Urick (6) formula. 
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Figure 2. The Dependent Relationship of a Upon V,~ at a Frequency of 1 MHz, 
and the Determination of the Value of the Constant P 


Key: (1) Experimental values 
(2) (In Urick's original graph, the vertical coordinates were 10 times 
larger) 
(3) Theoretical values 





Expression (23) shows that the sound phase velocity in the suspension is a 
function of frequency. Expression (24) shows that the sound absorption coef- 
ficient in the suspension is not diiectly proportional to the linear equation 
of frequency. These results can be explained by the Kramers-Kronig dispersion 
relationship (12) in quantum theory. 


In some natural coagulated suspensions, the dependent relationship of the sound 
phase velocity and the absorption coefficient upon frequency (for example, in 
clay sediments on the sea floor, C, and frequency are unrelated (13), a is 
directly proportional to the lineat equation of frequency) (13, 14) cannot be 
explained by the theoretical results described above. This shows that the re- 
sults of this article cannot be used for coagulated suspensions. The main 
mechanism of sound absorption in a coagulated suspension is internal friction 
of the suspension (10). 


VI. Comparison of Theory and Experiment 


Within the range of 0 «< V. < 0.4 of volume concentration of particles, Hampton 
(7) measured the dependent relationship of the sound phase velocity in the pure 
kaolinite suspension (frequency is between 50 kilohertz and 20 kilohertz) and 
the absorption coefficient (frequency is 100 kilohertz) upon the volume concen- 
tration of the particles. Within the frequency range of 10-10? kilohertz, Hamp- 
ton also measured the dependent relationship of the sound absorption coefficient 
in the pure kaolinite suspension upon frequency; the corresponding volume 








concentrations were 0.30 and 0.37, The results of the measurements are shown 

in Figures 3-5. The black dots in the diagrams are the experimental data of 
Hampton, the solid lines are the theoretical values of thie article, the dotted 
lines are the theoretical values of Ament, The percentages in Figure 5 are the 
volume concentrations of the particles, The physical parameters of the pure 
kaolinite syepension used were: 8) : 1,00 g/em*, pp = 2.71 /em’, r= 1,20 x 10cm, 
v= 0.01 em*/s, K; = 4.44 x 10°!? em*/dyn, Ky = 1.00 x 107!2 em® / ayn. 


Under Hampton's experimental conditions, equation (5) is satisfied. 
It can be seen from Figures 2-5 that: 


1. Consistency of the theoretical results of this article and the experimental 
data is very good. The theoretical results of this article are much better than 
the theoretical results of Ament. 


2. There is a visible frequency scattering effect in the pure kaolinite suspen- 
sion, and the velocity of sound increases as the frequency increases. 


3. The sound absorption coefficient in the pure kaolinite suspension is not a 
linear function of the volume concentration of the particles; when the volume 
concentration of particles is about equal to 0.17, the absorp.ion coefficient 
has a maximum value. 


4, Within a frequency range of 10 to 1,000 kilohertz, the sound absorption 
coefficient in the pure kaolinite suspension discussed is directly proportional 
to the nth degree equation of frequency; the degree n satisfies the inquality: 


1.14 <n < 1,82. 





The ratio of velocity in Kaolinite 
suspension to velocity in water 
versus volume concentration. 


XRX experimental data 
— calculated following Eq. (21) 
--- Calculated following Ament [3] 








Figure 3. The Dependent Relationship of C p/C) Upon V., a Comparison of 
Theory and Experiment 























”* 

















(pe 4 a —— = — 
10 * * v a a a a a eT 
VAM) frequency) kHe) 


Figure 4, The Relationship of ato V. Figure 5. Comparison of the Theoreti- 





at a Frequency of 100 kHz, cal and Experimental Rela- 
Comparison of Theory and tionship of a to Frequency 
Experiment ' 
Dependence of absorption coefficient Absorption coefficient as a function 
(in decibels perfoot) on volume con- of frequency. The volume concentra- 
centration at 100 kHz. tion of each kaolinite suspension is 
noted on each graph. 
. experimental data -.. experimental data 
— calculated following Eq (22) — calculated following Eq (22) 
— calculated following Ament — calculated following Ament 


VIL. Discussion 


In Reference (8) we only considered the mutual reaction of electrical forces 
between the particles in the jth shaded plane and the closest neighbors in the 
(j = 1)th and the (j + 1)th shaded planes, and we neglected the mutual reaction 
of electrical forces between that plane and the next closest and further neigh- 
bors. Actual calculations show that m,; and m) are numbers close to 1, and their 
contribution cannot be ignored. But because the value of P in Reference (8) was 
determined by experimental data, this shows that in Reference (8), m, and m» 
have already been unobtrusively added, i.e., the P determined in Reference (8) 
is numerically equal to H determined in this article. Although the value of P 
determined in Reference (8) was larger, equation (16) of Reference (8) is still 
effective. 


Experiments by Urick pointed out that the theory of this article can be applied 
within a volume concentration range of 0 < V. « 0.35 (see Figure 2), while the 
experimental results of Hampton indicated that the suitable range of the theory 
of this article can be even wider, i.e., it is applicable within the range of 
0<V_¢ 0.4. It seems that theoretical and experimental research in the theory 
of sound propagation in suspensions of a high concentration range is stili worth 
continuing. 
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LIFE SCLENCES 


ELECTRIC SHOCK THERAPY USED IN EPILEPSY TREATMENT 
Clinical Data 


Beijing ZHONGHUA SHENJING JINGSHENKE ZAZHI (CHINESE JOURNAL OF NEUROLOGY AND 
PSYCHIATRY] in Chinese No 1, 1981 pp 42-44 


[Article by Wang Zucheng [3769 4371 2110] and Yu Shicai [0060 0013 2088] of the 
Shanghai Psychiatric Hospital, Ma Huiying [7456 1979 3853] of Guizhou 304 Hospital, 
and Song Zhuanping [1345 1413 1627] of the Shanghai City Nanshi Ward Psychiatric 
Hospital: "Clinical Study of Electric Shock Therapy of Epileptic Mental Distur- 
bances; Analysis of 32 Cases") 


[Text] Electric shock therapy is a method of treating mental disturbance by arti- 
ficially inducing an epileptic reaction in order to better pinpoint symptoms of 
epileptic mental disturbance through electric shock treatment and thus effectively 
and rapidly control the symptoms and improve the effectiveness of treatment of men- 
tal disturbance accompanied by epileptic reaction. We have gathered for clinical 
analysis our hospital's information over the past 21 years on the use of electric 
shock to treat epileptic mental disturbances. 


lL. Source of Information 


Between August 1958 and June 1979, a total of 490 cases of epileptics were ad- 
mitted into our hospital. Electric shock treatment was used to treat epileptic 
mental disturbances in 32 of these cases, or 6.5 percent of the epileptics hospi- 
talized. In three cases, the patients were given repeated electric shock treat- 
ment when they were hospitalized again. Except in two cases, the nature of the 
epilepsy was recurring epilepsy. 


il. Analysis of Information 


|. General information: (1) There were 19 males and 13 females. Their ages 
ranged from 16 to 52 years old. In 7 cases the patients were 16 to 20 years old; 
in 13 cases the patients were 21 to 30 years old; in 8 cases the patients were 

31 to 40 years old; and in 4 cases the patients were over 41 years old. (2) The 
duration of hospitalization varied from 7 days to 14 years. There were 28 case- 
times in which hospitalization was less than half a year. (3) The history of the 
disease varied from 4 days to 20 years. In 17 case-times the case history was 
less than 5 years; in 8 case-times it was 5 to 10 years; in 8 case-times, 10 to 
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15 years; and in 2 case-times, 15 to 20 years. (4) In 18 case-times, mental 
disturbance occurred within 1 year; in 15 case-times, mental disturbance occurred 
within 1 to 5 years; in 2 case-times, mental disturbance occurred within 5 to 20 
years. In 18 of the 35 case-times, epileptic reactions occurred simultaneously 
with mental disturbance; for 16 inetances, epilepsy occurred before mental dis- 
turbance; and in 1 case mental disturbance occurred before epileptic reaction. 

(5) Clinical types of epilepsy and mental disturbances: There were 13 case-times 
(12 cases) of schizophrenic epileptic mental disturbance, 7 case-times (7 cases) 
of epileptic psychomotor reactions, 7 case-times (7 cases) of major and minor 
epilepsy accompanied by schizophrenia, 5 case-times (3 cases) of epileptic per- 
sonality disturbance, 2 case-times (2 cases) of mental disturbance caused by re- 
current epilepsy, and 1 case-time (1 case) of major epileptic reaction accompanied 
by menopause psychosis. (6) Combined medication: Except for 1 case, electric 
shock treatment was given in the remaining 34 case instances simultaneously with 
antipsychotic medication during hospitalization. The types of medical drugs in- 
cluded chlorpromazine, perphenazine, tardanum, fluoropiperidyl alcohol, reserpine, 
trifluorazine, fluoroperphenazine, and lithium carbonate. Daily dosage for treat- 
ment was within the commonly administered dosage. 


In 30 instances, electric shock therapy was combined with antiepileptic medication. 
Phenytoinum natricum was administered in 27 instances, phenobarbital was adminis- 
tered in 26 instances (in 21 instances, phenytoinum natricum and phenobarbital 
were administered in combination), primidonum was administered in 3 instances, 
meprobamatum was administered in 2 instances, and librium was administered in 1. 

In 5 instances, medication and shock therapy were not combined. 


2. Analysis of information on elect’ c shock treatment: (1) Electric voltage, 
electric current intensity, and duration of electric shock treatment: In the 
group of cases treated in combination with antiepileptic medicaiion, the threshold 
value of treatment of most cases was high, while in the group of cases treated 
without combination with antiepileptic medication, the threshold value was equal 
to the threshold value of ordinary treatment of psychosis. (2) The interval be- 
tween electric shock treatments and total number of treatments: In 16 instances, 
the interval of treatment was basically according to the usual interval of treat- 
ment--i.e., at the beginning, treatment was given two or three times a week, 
totaling five to eight times, and then treatment was changed to one or two times 
a week; 9 to 12 times of treatment were counted as one course of treatment. In 
another 15 instances, the treatment interval was often over 1 week, and there was 
no fixed pattern. In the other four cases, statistics were unavailable because 
the record of treatment was not detailed or the records were lost. In the total 
number of treatments, 6 treatments or fewer were given in 17 instances, 6 to 12 
treatments were given in 13 instances, and more than 12 treatments was given on 
only one occasion; there were 4 other cases with incomplete records. The total 
number of treatments was similar to the number given when treating schizophrenia. 
(3) The relationship betwee~ electric shock treatment and epileptic reaction: We 
marked out the month before electric shock treatment began and the month following 
the end of treatment. During this period, epileptic reactions occurred simultane- 
ously with the administration of electric shock treatment in 13 instances. During 
this period there were 22 instances in which there was electric shock treatment 
without its being accompanied by epileptic reaction. Among the 13 instances in 
which electric shock treatment was given at the same time as the occurrence of 
epileptic reactions, the treatment stimulated an increase in epileptic reactions 
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in 4 instances, and the treatment produced a decrease in epileptic reactions in 
9 instances. 


3. Analysis of the effects of electric shock therapy: (1) The effectiveness of 
treatment of various clinical types: Electric shock therapy was effective in 5 of 
? instances of psychomotor epileptic reaction, in 9 of the 13 instances of schizo- 
phrenic epileptic psychosis, in 4 of the 7 instances of major or minor epileptic 
reactions accompanied by schizophrenia, and in 2 of the 5 instances of epileptic 
personality disturbance. It was effective in one of the three instances of other 
epileptic mental disturbance or epileptic reactions accompanied by psychosis. Of 
the 35 case-times, it was effective in 21 instances, and so the total percentage 
of effectiveness was 60 percent. It was ineffective in both of the cases of recur- 
rent epilepsy. (2) The effectiveness of treatment of individual psychotic symptoms: 
Electric shock therapy was effective in the treatment of symptoms of speaking to 
oneself, disorderly behavior, melancholia, disorganized thought, excitability and 
being easily angered, impulsive destruction, and delusion, as shown in the accom- 
panying table. (3) The relationship between effectiveness of treatment and the 
course of epilepsy and psychosis: There were 18 instances of simultaneous occur- 
rence of epilepsy and mental disturbance; treatment was effective in 13 of these, 
or 72.2 percent. There were 16 cases of epilepsy before the occurrence of mental 
disturbance; treatment was effective in 8, or 50.0 percent. There was one in- 
stance in which mental disturbance was followed by epileptic reaction, and treat- 
ment was ineffective. It seems that the effectiveness of treatment of the first 
type of cases was becter than in the second type of cases, while there was only 
one case in the third type and judgment could not be made. (4) Relationship be- 
tween effectiveness of treatment and the frequency of treatment: In 9 instances, 
effectiveness of the treatment occurred within 6 treatments, while in 8 instances 
the treatment was ineffective; thus the percentage of effectiveness was 52.9 per- 
cent (9/17 case-times). In 10 instances, effectiveness occurred after 7 to 12 
treatments, while in 3 instances the treatment was ineffective; thus the percen- 
tage of effectiveness was 76.9 percent (10/13 case-times). The percentage of ef- 
fectiveness is higher [when the frequency of treatment is greater] than otherwise. 


Accompanying Table. Effectiveness of Electric Shock Therapy in the Treatment of 
Psychotic Symptoms 
Key: 
A. Psychotic symptoms 
B. Case-times 


(A) (B)(C) (D) (A) (B)(C) (D) C. Effective case-times 
— * D. Percentage of effectiveness 
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4. Side effects after electric shock therapy: In 35 case-times, following some 
250 electric shock treatments, there were 6 cases (7 case-times) that showed side 
effects, or 20 percent of the total case-times. Side effects were manifested as 
follows: In four cases, major epileptic reactions occurred frequently after two 
electric shock treatments. In one case, semiconsciousness, excitation, and a manic 
state accompanied by refusal to eat occurred after six electric shock treatments; 
these effects persisted for 4 days, and then a major epileptic ‘ceaction »ccurred 
and the condition improved. In 1 case, a seubstiffness state occurred after 10 
electric shock treatments and the condition improved later by itself. 


IIl. Discussion 


There have already been reports on the use of tic to treat epileptic reactions. 
The reaction threshold of epilepsy can be increaced in epileptic patients who have 
used certain medicines (such as xetrazol) to cause convulsions. The reports have 
also mentioned frequent death due to recurrent epilepsy of "combined epileptic and 
psychotic symptoms.” If this type of reaction can be induced by man artificially, 
the duration of mental disturbance of the patients can be greatly shortened. 
There...re, Kalinowsky believed: Artificially induced epilepsy serves a protective 
function in automatic epilepsy. If antiepileptic drugs still cannot control epi- 
leptic reactions, electric shock can be used to supprese such reactions. He also 
pointed out: Incomplete epileptic reactions can .e terminated by one or two 
automatic or artificially induced epileptic reactions. Gerletti mentioned that 

the use of electric shock therapy to treat epileptic patients can suppress auto- 
matic epileptic reactions for a longer period of time. Some writings have also 
pointed out that the use of electric shock therapy to treat epileptic mental dis- 
turbance, psychomotor reactions, schizophrenic epileptic menial disturbance, and 
somnolent states has a definite effect, and it can check some majo~ epileptic 
reactions and emotional disturbance and improve the condition. Ouchi et al be- 
lieved: "Electric shock treatment of schizophrenic epileptic delusion symptoms 

4 to 10 times has a definite effect, and it can also eliminate hallucinations, 

but there are also ineffective cases." 





From our indications, the most suitable types of epileptic mental disturbance for 
electric shock treatwent are: psychomotor reactions, personality disturbance, 
schizophrenic mental disturbance, and major and minor epilepsy accomp. nied by 
schizophrenic patients. The effectiveness is always 60 percent and close to the 
percentage of effectiveness (66 percent) in the treatment of schizophrenia by 
electric shock therapy mentioned by Cornelio. In 13 cases, major epilepsy oc- 
curred within th. period from 1 month before to 1 month after electric shock 
treatment. Of these, there were nine cases (69.2 percent) in which electric 
shock treatment caused a reduction in epileptic reactions; this is consistent 
with reports that electric shock therapy can suppress major epileptic reactions. 
But we should also note that in a few cases (four cases, 30.8 percent), treatment 
was terminated because of more frequent reactions. In addition, we also note a 
total of six ceses of transient semiconsciousness, excitation and manic activity, 
and substiffness state wnich subsided and the condition of the patients improved 
without any visible after effects. Thus, electric shock therapy is a relatively 
safe method of treating automatic epileptic mental disturbance. 
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in Je cases (35 Casewtimes), 30 were given antie,s.leptic medication simultaneously 
ind Che treatment threshold values were mostly higher. The treatment threshold 

lue in the five instances in which medication was not given was equal to the 
thveshold values of electric shock treatment of other kinds of paychosis. There- 
‘ore, epllepsy does not increase the treatment threshold value, but the threshold 
values increased only after taking antieptieptic medicines. 


in the method of electric shock therapy for treatment of epileptic mental dis- 

irbance, the interval of treatment and the total frequency of treatment are basi- 
cally the same as those in ordinary electric shock therapy. Although nearly one- 
valf (48.5 percent) of the cases in the data received a total of six or fewer 
(reatments, the effectiveness rate was only 52.9 percent. The effectiveness of 
(reatment was 76.9 percent for cases treated 7 to 12 times, far higher than those 
cases treated six or fewer times. It can be seen that to obtain the maximum re- 
ult of Creatment, the total number of customary electric shock treatments is re- 
quired, 


ibe mechaniem of electric shoe therapy for treating epileptic mental disturbance 
has been explained by Kalinowsky: Electric shock treatment can temporarily in- 
crease the tic threshold value. From this we can infer that electric shock in- 

reases the discharge thr shold value manifested as psychotic symptoms and thus 
-auses the psychotic symptoms to disappear. Because we still are unable to ex- 
plain epilepsy and the mental disturbance caused by it perfectly, and because the 
mechaniem of electric shock therapy is not completely clear, we can only make such 
an inference. 


Clinical Analysis 


Seijing ZHONGHUA SHENJING JINCSHENKE ZAZHI [CHINESE JOURNAL OF NEUROLOGY AND 
PSYCHIATRY] in Chinese No 1, 1981 pp 45-46 


Article by Lu Longguang [7627 769% 0342], Ma Wenhe [7456 2429 0735], Zhane Xinbao 
1728 1800 0202), and Zhao Kiaoli [6392 2556 7787) of the Nanjing Psychiatric 
fospital: “The Experience of Outpatient Electric Shock Treatment; Clinical Anal- 
sis of 800 Cases") 


fext] Electric shock therapy already has a history of 40 years. Its active 
funetion in psychiat fe treatment has been affirmed by the world's nations for a 
ong time. To expiore the clinical value of electric shock therapy before admin- 
istering antipsychotic medication, we reviewed 1,542 cases of various types of 
osychosis which received outpatient electric shock therapy within the 3 years 
from 1952 to 1954, Of these, information on 800 cases was complece. ihe follow- 
ing is a report of the effectiveness of treatment and a review of some of the 
ases: 


Analysis of Information 


Sex and age: In tnis group, 412 cases were males, and 388 cases were femaies. 
ight cases were less than 15 years old, 91 cases were between 16 and 20 years 
‘id, 389 cases were between 21 and 30 yeare old, 178 cases were between 31 and 40 
years old, 93 cases were between 41 and 50 years old, 36 cases were bt tween 51 and 
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60 years old, 2 cases were between 6] and 70 years old, and 1 case was ov © 70 
years old, The ages of the patients receiving treatment were between 13 and 71 
years old; 650 cases ($1 percent) were between 21 and 50 years old. Those under 
15 and over 50 totaled 49 cases. 


2. Prequency of treatment: This group received a total of 6,543 electric shock 
treatments. Each person received an average of 8.1 treatments; the maximum number 
of treatments given to one patient was 38; the minimum number of treatments given 
to one patient was one, 


3. Effectiveness of treatment: 


(1) Type of disease and effectiveness of treatment within a short period: In the 
various types of psychosis, the best result was in the treatment of melancholia 
and mania, the percentages of cure in these diseases being 89.9 percent and 86.7 
percent, respectively. The percentage of cure for affective psychosis was 82.2 
percent. These were for the most part acute cases. While they were being admin- 
istered electric shock .herapy, moet cases received psychotherapy. The percentage 
of cure for compulsive neurosis was 44.4 percent, while for schizophrenia it was 
40.4 percent (Table 1). 


Table 1. Type ot Disease and Effectiveness of Treatment Within a Short Period 
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Key: 

Types 

Schizophrenia 

Mania 

Depression 
Affective psychosis 
Compulsive neurosis 
Total 

Cured 

Number of cases 
Visible improvement 
laprovement 
Ineffective 

Total 


— 
owoeoxw Ow & wwe 


— = bee 
why —— 
—- = . 





(2) Duration ef divease and effectiveness of treatment within a short period: The 
percentage of patients cured within 1 year of the disease was 67.2 percent, the 
percentage of patients cured within | to 2 years was 59.9 percent, the percentage 
of patients cured within 2 to 3 years wae 41.8 percent, the percentage of patients 
cured within 3 to 5 years was M& percent, and the percentage of patients cured 
after 5 years of the disease was 32 percent. The shorter the history of the 
disease, the higher the effectiveness; this (finding) is consistent with domestic 
and foreign reports, The percentage of cure for patiente with chronic disease of 
over 2 years’ duration is stil! optimistic (Table 2). 


lable 2, Duration of Disease and Effectiveness of Treatment Within a Short Period 
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Key! 
1. Duration of disease (year) 
2. Cured 
3. Number of cases 
4. Visible improvement 
5. Improvement 
6. Ineffective 
7. Total 


4, Complications: Among the complications during the period of treatment of this 
group, failing memory was the most commonly seen, there being 28 cases or 3.5 per- 
cent. This complication occurred mostly after three to eight treatments. Head- 
aches occurred in 25 cases (3.1 percent), all during the first few treaiments. 
Spinal pain occurred in 23 cases (2.9 percent). Physical and X-ray examinations 
revealed one case (0.1 percent) of compressed fracture of the thoracic vertebrae. 
Bleeding of the gums occurred in 22 cases (2.8 percent). Injury of bitten tongue 
oceurred in 14 cases (1.8 percent). Pain in the shoulder and knee joints occurred 
in 14 cases (1.8 percent). Mental disturbance occurred in 8 cases (1 percent) 
(including semiconsciousness, acute manic activity, and anxiety which were psy- 
chotic symptoms not seen in the original disease). Breathing difficulty occurred 
in 5 cases (0.6 percent), loosening of the front teeth occurred in 2 cases (0.3 
percent), and aspiration pneumonia occurred in 1 case (0.1 percent). 


The number of complications as reported was not uniform. The causes of the com- 
plications of this group included not only the effects of the treatment itself, 
but also a lack of technical skill, insufficient pretreatment preparation, and 
failure to follow normal procedures. 








The question of deaths: Since our hospital's outpatient service began using elec- 
tric shock therapy, several tense of thousands of treatments have been admin.stered, 
and there has not been a single case of death. Thus it can be seen that this is 
still a relatively safe type of physical therapy. 


>. Long-term effectiveness of electric shock therapy for schizophrenia: Of the 
526 cases of schizophrenia, 190 cases reported results of followup visits accord- 
ing to the outpatient followup visit and observation records. in some of the 
cases, aftereffects of insomnia and such symptoms of neurosis remained after con- 
ditions improved. Except for a few patients of this group who regularly took 
sedatives and sleeping pille, no other treatment was given. Of the 53 cases of 
followup visite within 1 year, 67.4 percent remained cured or showed visible im- 
provement. Of the 44 cases of followup visite 3 to 5 years after treatment, 47.7 
percent remained cured or nearly cured (Table 3). 


Table 3, Effectiveness of Followup Treatment of 190 Schizophrenic Cases 


























> 
(6)% tt 9 jie | 

Key: 

1. Time of followup visit (year) 

2. Cured 

3. Visible improvement 

4. improvement 

5. Ineffective 

6. Total 

7. Number of cases 
Discussion 


Electric shock is a kind of physical therapy. Because the modus operandi of this 
kind of treatment is special, and because of the strong reaction of the patients, 
evaluation [of the treatment) has not been consistent since its inception. But 
because of its true effectiveness, it has thus been used continuously. This 
therapy was once banned from use during the latter part of the 1960's and the 
early 1970's in our nation. At present, most hospitals have revived its use, but 
the number of cases treated and the frequency of application have lessened. The 
same situation exists abroad. To explore the clinical value of electric shock 
therapy, we have summarized and reviewed past work and have proposed some views. 
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|. A@ the methods of treatment of psychosis continue to improve, the effective- 
jesse of treatment has aleo gradually increased, especially during the past more 
than <U years, The emergence and the rapid development of antipsychotic medicines 
have opened up a wide tield for the treatment of paychosis, replacing most of the 
functions of eleetric shock therapy. Overall measurement of the information of 
this group shows that the method of electric shock is etill an effective treatment 
for some types of psychosis, especially for melancholia and serious suicidal cases, 
Au.tipeyehotic drugs can be used in coordination with it, but they cannot replace 
ite function, 





2. Of the 800 cases of paychotic patients subjected to electric shock therapy, 
54.8 percent were cured by short-term treatment, and the effectiveness was 83 per- 
cent. The effectiveness of treatment of melancholia wae 95.5 percent, the effec- 
tiveness of treatment of mania was 94.7 percent, and the effectiveness of treat- 
ment of schizophrenia was 76.9 percent. Outpatient cases followed up for between 
half a year and 5 years numbered 190; the percentage of those cured was 67.4 per- 
cent, and the effectiveness was 79 percent. This figure is definitely convincing. 
Foreign reports have pointed out that both before and after antimelancholia drugs 
came on the market, electric shock therapy has always been an outstandingly effec- 
tive method for ameliorating melancholia, controlling suicidal tendencies, and 
consolidating the effectiveness of treatment. Facing the threat of serious sui- 
cidal tendencies, it can rapidly improve and eliminate the symptoms. Some people 
have dubbed it a “Lifesaving treatment." Hoffet pointed out that the effective- 
ness of treatment of melancholia by promethazine is 54.5 percent, while the effec- 
tiveness of electric shock therapy is 95 percent. The effectiveness of treating 
schizophrenia with a touch of melancholia is 83.5 percent, which is a higher per- 
centage than if only antipsychotic drugs are used. The complications of electric 
shock therapy are no greater than those in using antipsychotic drugs, and the 
death rate is the same. Negating electric shock therapy too early lacks suffi- 
cient factual basis. 


We believe that continued use of electric shock therapy should require a strict 
selection of the cases, so as to grasp the suitable indications and contraindica- 
tions. A detailed and careful physical examination should be carried out prior 
to treatment in order to profoundly grasp the overall situation of the advantages 
and the shortcomings. Patients with serious suicidal tendencies should be consid- 
ered for electric shock therapy first if there are no contraindications. Refusal 
to eat, stiffness, psychomotor excitation with symptoms of aggressive attack on 
other people and self-mutilation which are difficult to control, and catatonic 
schizophrenia accompanied by a touch of melancholia are all suitable indications 
for this method of treatment. Information of this nature shows that patients who 
react well to electric shock therapy generally will respond well to treatment 
after two or three treatments. If there is still no improvement after more than 
five treatments, it will be difficult to obtain satisfactory results if treatment 
is continued. Also, some complications such as mental disturbance and memory 
disturbance all occur after eight treatments; this coincides with the reports in 
medical documents concerning the direct proportion between the seriousness of 
memory disturbance and the frequency of electric shock treatments. The patient 
being older than 50 years of age, relatively poor physical health, or psychosis 
accompanied by cardiovascular disease should be contraindications. Young people 
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under 15 years of age have not fully matured, and such treatment should be strictly 
limited, According to reports, analysis of the cause of death in electric shock 
therapy reveale that cardiovascular disease, prostrated breathing, and cerebro- 


vascular disease are the three most frequently seen causes. These can serve as 
reterences, 
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